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Abstract. We numerically study the orientation deformations in nematic liquid crystals around charged 
particles. We set up a Ginzburg-Landau theory with inhomogeneous electric field. If the dielectric 
anisotropy ei is positive, Saturn ring defects are formed around the particles. For ei < 0, novel "ansa" 
defects appear, which are disclination lines with their ends on the particle surface. We find unique defect 
structures around two charged particles. To lower the free energy, oppositely charged particle pairs tend 
to be aligned in the parallel direction for ei > and in the perpendicular plane for ei < with respect 
to the background director . For identically charged pairs the preferred directions for ei > and £i < 
are exchanged. We also examie competition between the charge-induced anchoring and the short-range an- 
choring. If the short-range anchoring is sufficiently strong, it can be effective in the vicinity of the surface, 
while the director orientation is governed by the long-range electrostatic interaction far from the surface. 

PACS. 61.30.Dk Continuum models and theories of liquid crystal structure - 61.30.Jf Defects in liquid 
crystals - 77.84.Nh Liquids, emulsions, and suspensions; liquid crystals - 61.30.Gd Orientational order of 
liquid crystals; electric and magnetic field effects on order 
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1 Introduction 

A variety of mesoscopic structures have been found in liq- 
uid crystals around inclusions such as colloids and water 
droplets In nematics, inclusions distort the ori- 

entation order over long distances, inducing topological 
defects [J.FJT.'alElliniiniinillS]. We mention the forma- 
tion of structures or phases, such as string- like aggregates 
[2l lT2llT4l[T5] . soft sohds supported by a jammed cellular 
network of particles |16j , and a transparent phase includ- 
ing microemulsions |17lll8j . The origin of the long-range 
distortions has been ascribed to the anchoring of the liq- 
uid crystal molecules on the inclusion surface [HISlinjEllHlini 
[Tgi[TH[ni[T5|[Tg|[2U]. It arises from the short-range molec- 
ular interactions between the liquid crystal molecules and 
the surface molecules. In the Ginzburg-Landau-de Gennes 
theory, we have a surface free energy depending on the ori- 
entation of liquid crystal molecules on the surface. 

In this paper, we are interested in another anchoring 
mechanism. That is, electrically charged inclusions align 
the liquid crystal molecules in their vicinity to lower the 
electrostatic energy [^[^[^ . which can be relevant for 
ions and charged particles. In fact, de Gennes at- 
tributed the origin of the small size of the ion mobility 
in nematics to a long-range deformation of the orienta- 
tion order around ions. However, the effect of charges in 
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liquid crystals remains complicated and has rarely been 
studied, despite its obvious fundamental and technologi- 
cal importance. It is of great interest how the electric- field 
anchoring mechanism works and how it is different from 
the usual short-range anchoring mechanism. 

The electric field and the liquid crystal orientation are 
coupled because the dielectric tensor depends on the 
local orientation tensor Qij (see equation (14)). The align- 
ment along a homogeneous electric field is well-known |25j , 
but the alignment in an inhomogeneous electric field has 
not yet been well studied. When the dielectric tensor is 
inhomegeneous, it is a difficult task to solve the Pois- 
son equation and seek the electric potential 'P. We here 
perform numerical simulations placing charged particles 
in liquid crystals in a three-dimensional cell. We use the 
Ginzburg-Landau-de Gennes scheme in terms of the ori- 
entation tensor Qij [23,25,26,27J. A similar approach has 
recently been used to calculate the polarization and com- 
position deformations around charged particles in elec- 
trolytes It is worth noting that hydration of water 
molecules around ions is analogous to the orientation an- 
choring of liquid crystal molecules around charged parti- 
cles, as pointed out by de Gennes [24j[25|. 

In Section 2, we will present a Ginzburg-Landau-de 
Gennes theory for liquid crystals containing charged par- 
ticles. In particular, we will give two general forms of the 
electrostatic free energy for the fixed-charge and fixed- 
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potential cases (which can be used for any dielectric fluids 
containg charges). In Section 3, we will explain the nu- 
merical method adopted in this work, in Section 4, we will 
present numerical results of equilibrium configurations of 
the orientation order around charged particles. We will 
also examine competition of the short-range and electric- 
field anchoring mechanisms. In Section 5, a summary and 
critical remarks will be given. 



Charge Qt Potential * ^ 
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{a) Change fixed 
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(b) Potential fixed 



2 Theoretical Background 

We consider a liquid crystal system in a cubic box and 
place one or two charged spherical particles with radius R 
inside the box. The particle positions are written as i?„ 
(n = 1, 2). The liquid crystal order is described in terms of 
the symmetric orientational tensor Qij(r) with the trace- 
less condition Qu = [53]. We place one or two charged 
particles with radius R considerably longer than the ra- 
dius of the solvent molecules. In this work the Boltzmann 
constant is set equal to unity and the temperature T rep- 
resents the thermal energy of a liquid crystal molecule. 



2.1 Model 

We are interested in the equilibrium liquid crystal orienta- 
tion around the particles, which minimizes the sum of the 
Landau-de Gennes free energy, the short-range anchoring 
energy, and the electrostatic energy. Thus the total free 
energy of the liquid crystal containing charged particles 
consists of four parts as [2l ][22l[23] 



J- — J-Q + J-g + J-a + J^e- 

The first term is of the Landau-de Gennes form, 



^0 



dr 



where we introduce 



J2 



)2. 



^ I + ^ 



J3 — QijQ jkQki 



(1) 



(2) 



(3) 



Hereafter repeated indices are implicitly summed over. 
The coefficient A is dependent on the temperature T, 
while the coefficients B and C are positive constants as- 
sumed to be independent of T. The second term is the 
gradient free energy in the one-constant approximation. 



(4) 



where Vfc — d/dx^ {xk — x, y, z) are the space derivatives 
and L is a positive constant. The space integrals / dr in 
equations (2) and (4) are to be performed only outside 
the particles \r — > R. It is convenient to define the 
length, 

d = T/L, (5) 



Fig. 1. A capacitor and an inhomogeneous fluid containing a 
net charge Qin in the fixed-charge case (a) and in the fixed- 
potential case (b). The charge and potential of the lower plate 
are Qb and <5b, while those of the upper plate are Qt and $t- 



which is the typical molecular size of liquid crystal. The 
term Ta represents the short-range anchoring free energy. 
It is expressed as the integral on the particle surfaces. 
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where da is the surface element, f is the outward normal 
unit vector to the surface, and w represents the strength 
of the anchoring. For the uniaxial form Qij = S{ninj — 
5ij/i), we have Ta — wS J da[l/3 — {v ■ n)^]. Thus, for 
neutral particles, positive and negative values of w lead to 
homeotropic and planar anchoring, respectively. 

We explain the electrostatic part Te, which depends 
on the experimental method. As a generalization of the 
theory by one of the present authors [2T], we allow that 
the fluid region can contain a net charge Qin = J drp, 
where p = p(r) is the charge density inside the fluid. As 
in Figure 1, we insert the fluid between parallel metallic 
plates in the region Q < z < H . The surface charge and 
the potential of the lower plate at z = are Qb and ^b, 
while those of the upper plate aX z = H are Qt and 
We require the overall charge neutrality condition. 



Qin + Qb + Qt = 0, 



(7) 



since the eletric field in the metal plates should vanish. In 
terms of Q = (Qt — Qb)/2, we may set 



Qt 



-Q - Qin/2, Qt = Q - Qin/2. 



(8) 



(i) In (a) in Figure 1, Q can be fixed and can be a control 
parameter, where the potential difference. 



(9) 



depends on the fluid inhomogeneity induced by the chaged 
particles. Here the electrostatic energy of the surface charges 
of the plates is fixed, the appropriate form of is 



1 

8^ 



drD E 
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where Ei = —d^l'/dxi is the electric field and Di — SijEj 
is the electric induction with Sij being the dielecric tensor. 
Here we superimpose small variations (5Qb, ^Qt, ^P, and 
5£ij on Qb, Qt, P, and e^, respectively. We use the relation 
/ drE ■ SD/Att = -PbSQb + <PtSQt 
obtain the incremental change of !Fe 



J dr'PSp. We then 



as 



5Te = VSQ+ / drSp[<P 



1 

'8^ 



drdSijEiEj. 



(11) 



(it) On the other hand, in (b) in Figure 1, the potential 
difference V can be fixed and can be a control parameter 
with Q being dependent on the fluid inhomogeneity. The 
appropriate form of J-'e is 



1 

8^ 
= f dr 



drD -E-VQ 



pl-P- 



D E 



(12) 



where the second line follows from the second line of equa- 
tion (10). This is the Legendre transformation of the elec- 
trostatic free energy in the fixed-charge case. Here we use 
the same notation J-'e in the two cases. Then the incre- 
mental change of !Fe reads 



-QSV^ 
1 

'8^ 



drSpi <P 



drSsijEiEj. 



(13) 



where the first term on the right hand side is different from 
that in equation (11). It is worth noting that the second 
line of equation (12) yields the frequently used expression 
J^e = — J drD ■ E/8tt in the fixed-potential condition for 
dielectric fluids without charge {p — 0) (see reference [24. . 
for example). 

The potential 'P satisfies the Poisson equation, 

V,(eyV,<Z') = -Anp. (14) 

We assume the linear form of the dielectric tensor, 



eij{r) = EoSij + eiQij{r), 



(15) 



in the liquid crystal region (the particle exterior) Defin- 
ing in the whole space, we may solve equation (14) by 
setting Sij{r) = £pSij in the particle interior. Then the in- 
tegrals in equations (10) and (12) are over the whole cell 



region. We then have SJ^e/SQij 
fixed Q and at fixed V. 



— —EiEiEj/Sn both at 



^ In the nematic state we have ey = eo + 25'ei/3 along the 
director n and = eo — Sei/3 in the perpendicular directions 
[24j . where the amplitude S is given in equation (20). 



2.2 Equilibrium conditions 

In our numerical work we will adopt the geometry (b) in 
Figure 1 and set = 0. The charge density p is fixed. We 
define the tensor, hij = dJ-/SQij + XSij, where A is chosen 
such that hij becomes traceless. Some calculations give 



"ij — {A + CJ2)Qij — B ( QikQkj — 2>^2% 



£l 



EiEj 



In equilibrium, minimization of J- yields 



hi 



0, 



(16) 



(17) 



in the particle exterior. The boundary condition of Qij on 
the particle surface is given by 



Lf ■ S/Qij + w{viVj — Sij/i) — 0. 



(18) 



Obviously, the defect structure is independent of the sign 
of the particle charge, since Qij is coupled to the bilinear 
terms of E in equation (16). 

For B > Q uniaxial states are selected in the bulk re- 
gion below the isotropic-nematic transition A < At [5S], 
where Qij — S^niUj — Sij/S) and 



At = B^/27C. 



(19) 



Substituting the uniaxial form into the first line of equa- 
tion (16), we obtain 2CS'^-BS + 3A = 0, which is solved 
to give 



B/4C + [{B/ACf - 3A/2C]i/^ 



(20) 



Just below the transition we have S—St = B jZC . How- 
ever, it is known that the liquid crystal order is consider- 
ably biaxial inside defect cores [51 [5^l26j . See Figure 3 of 
Ref.[23j for the biaxiality of the Saturn ring core (where 
the spatial mesh size is finer than in this work) . Note that 
Qij can generally be expressed as 



Si{ninj - 5ij/i) + S2{mim-j - tilj), 



(21) 



where n, m, and £ constitute three orthogonal unit vec- 
tors. Inside defect cores, the amplitude 5*2 of biaxial order 
is of the same order as the amplitude 5*1 (= S in this work) 
of uniaxial order. Outside the defect cores, 5*2 nearly van- 
ishes and the orientation order becomes uniaxial. 

In addition, the polarization vector of the liquid crystal 
is given by Pi = Xij Bj in terms of the susceptibility tensor 
All Pi, we have 



Xij- From EijE^ = E, 



(22) 



This tensor should be positive-definite in equilibrium to 
ensure the thermodynamic stability in the (paraelectric) 
nematic phase For the special form (15) this require- 
ment becomes 

eo - 1 + eiqa > 0, (23) 
where qa (a = 1, 2, 3) are the eigenvalues of Qij. 
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2.3 Electric field effect near the surface 

Let us consider the electric field effect near a particle sur- 
face. For simplicity we assume |ei| <£q- Then the sur- 
face electric field Es is estimated to be of order eZ/eoR'^, 
where Ze is the particle charge (with e being the elemen- 
tary charge), (i) Far above the transition A ^ At in the 
isotropic phase, we neglect the terms proportional to B, 
C, and L in equation (17) to obtain Qij = eiE'^{xiXj /r"^ — 
5ij /Z) /?>nA, which grows as A is decreased as a pretransi- 
tional effect, (ii) Just below the transition, a nonlinear de- 
formation occurs for |ei|i?g/87r >AtSt — B'^ /%\C'^, which 
is easily realized for small B. (iii) In the nematic phase 
far below the transition, strong nonlinear deformations of 
Qij are induced on the surface for R < £ with [23] 

£^\Z\{\ei\eBd/127reoSy/^ (24) 

where d is defined by equation (5) and 

£b = eyeoT (25) 

is the Bjerrum length. This criterion arises from the bal- 
ance of the gradient term (~ LSR^^) and the electrostatic 
term (~ EiE'^/Stt oc i?^^) in hij in equation (17). Further- 
more, for sufficiently large (/R, a defect is formed around 
the particle, where the distance from the surface is of order 
i-R. 

It is important to clarify the condition of defect for- 
mation in real systems. Let us assume Eq ^ 2, |£i| ^ Eq, 
5* ~ 1, d ~ 2nm, and ig ^ 24nm. Then £ ^ |Z|nm. Thus, 
the relation R < i holds for microscopic ions, though our 
coarse-grained model is inaccurate on the angstrom scale. 
See the remark (3) in the last section for a comment on 
ions in liquid crystal. We may also consider a large particle 
with a constant surface charge density 

(7 = Z/4:TtR^. (26) 

It may be difficult to induce sufficient ionization on col- 
loidal surfaces in liquid crystal solvents. One method of 
realizing charged surfaces will be to attach ionic surfac- 
tant molecules on colloidal surfaces. For such a particle, 
the condition of defect formation becomes R^ Rc , where 

i?, = {3eQS/47T\ei\£Bdy/^cr~\ (27) 

Using the above parameter values, we have Rc O.lcr^^nm 
(with (7 in units of nm^^). For example, if cr = 0.0624nm^^ 
or ea — 1/iC/cm^, we obtain Rc — 1.6nm. Here the elec- 
tric field at the surface is ea/iweo ~ 100V//im, which 
is strong enough to align the director field. Electric field 
applied macroscopically is typically of order lV//j,m [501 

IS]. 

3 Simulation method 

We give our simulation method in the Landau-de Gennes 
scheme under the condition of V ~ 0. For simplicity, we 







1 ' 1 

dJ2/dr ~ 

/g 







1 2 

ir-R)/R 



Fig. 2. Derivative dJ'^/dr in units of dT^ and gradient free 
energy density fg = L{\7kQij)^ /2 in units of Td~^ vs normal- 
ized distance (r — R) / R from the surface of a charged spherical 
particle. The path starts from a surface position and passes 
through a Saturn ring (see Figure 3). 

impose the periodic boundary condition in the xy plane. 
We suppose nanoscale particles confined between a thin 
layer. 

In the previous section we have assumed sharp bound- 
aries between the particles and the liquid crystal region. 
However, precise simulations are not easy in the presence 
of sharp curved boundaries on a cubic lattice, unless the 
mesh size is very small. In this work, to overcome this dif- 
ficulty, we employed the smooth particle method. That is, 
we introduce diffusive particle profiles by [Sl fTTlfT^iroiEn] 

0n(r) = ^tanh{ ^-'^"-' } + i, (28) 

where the surface is treated to be diffuse with thickness 
d = T/ L in equation (5), i?„ represents the particle cen- 
ter, and R is the particle radius. 

In terms of (t>n{T)j the overall particle and charge dis- 
tributions are expressed as 

0(r) = ^</.„(r), (29) 

n 

p(r)= -^Z„<^„(r), (30) 

n 

where Z„e are the particle charges and v — ^.-kR? ji is 
the particle volume. The charge distribution is assumed 
to be homogeneous inside the particles. In Tq in equation 
(2) and Tg in equation (4), the space integrals outside the 
particles j' dr should be redefined as 

/ rfr-(...) = 1 dr[l-0(r)](..-). (31) 
The surface integral in equation (6) is also redefined as 

y da(...)= y dr|V0|(---). (32) 

Then, the short-range anchoring free energy (6) is rewrit- 
ten as 

Ta = -w I drg,,(V,(/.)(Vj(/.)/|V0|. (33) 
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The dielectric tensor is space-dependent as 
eij{r) = [eo + (sp - £Q)4>\5tj + ei(l - 



IJ7 



(34) 



where Ep is the dielectric constant inside the particles. 

To seek Qij satisfying equations (17) and (18), we 
treated Qij{r, t) as a, time-dependent tensor variable obey- 
ing the evolution equation, 




ld_ 



(r,t) = - 



-(1-, 
w 



)/t,j -i(Vfe0)(VfeQ.y) 



1 2 ^ij 



Fig. 3. (color online) Orientational field around a charged 
particle for (a) Z = 60, (b) Z = 100 and (c) Z = 160 in a 
nematic solvent with ei — l.Seo- Short lines (in blue) represent 
the director n = (n^, ny,n^) and cylinders (in green) in (b) and 
(35) (c) contain a Saturn ring. 



where C is a constant kinetic coefficient. In the first line, 
the functional derivative is taken both inside and outside 
the particles with the redefinitions (29)-(34), with A ensur- 
ing Qn = 0. In the second line, hij is defined in equation 
(16) and Vk4> arises from the factor 1 — </> in equation (31). 
On a cubic 64 x 64 x 64 lattice, we integrated the above 
equation for Qij . Space and time are measured in units of 
d and 

r = (f/QL, (36) 

respectively. The space mesh size is d and the time mesh 
size is At — O.Olr in the integration. The cell interior is in 
the region ^ x, y, z ^ 64(i. We solved the Poisson equa- 
tion (14) at each integration step using a Crank-Nicolson 
method [25] . 

As the boundary conditions of Qij at z = and QAd, 
we assume the homeotropic anchoring Ui = Siz for ei > 
and the parallel alignment tii = 5ix for Si < 0, where 
n = (nx,ny,nz) is the director with i — x,y,z. Those of 
Qij in the x and y directions are the periodic boundary 
conditions. The potential <P vanishes at z = and 6Ad 
and is periodic in the xy plane. Note that the electric field 
at z = and GAd is along the z axis, so the electrostatic 
energy is lowest for the selected director alignments both 
for £i > and £i < 0. In order to approach a steady 
state, we performed the integration until \dJ-/dt\ became 
less than IQ-^T/t. 

In our steady states thus attained, we confirmed that 
both equations (17) and (18) excellently hold in the bulk 
liquid crystal region and near the particle surfaces, re- 
spectively. Mathematically, they should hold in the thin- 
interface limit d <C i?, where — = 6{r — R)v around 
a spherical surface with u being the normal unit vector. 
In Figure 2, we show our numerical result of the deriva- 
tive dJ^jdr = 2QijdQij/dr and the gradient free energy 
density fg = L{V kQijY /"^ around a particle surrounded 
by a Saturn ring defect for w = 0. See the next section 
for details of the calculation and Figure 3 for its 3D pic- 
ture. We can see that dJijdr is nearly equal to zero at 
the surface and exhibits double peaks around the Saturn 
ring position. The boundary condition v ■ '^Qij = in 
equation (14) is thus nearly satisfied even in the presence 
of a defect in our diffuse interface model. 



4 Numerical results 

In our simulations, we set 

v4=-15T/d3, B= 1^1/2, C = 3B, 
Ib = 12d, = eo- 

For example, for Eq = 2.3 and T ~ 300K, we have d = 
2nm, Ib = 24nm, and L = 2pN. The nematic order pa- 
rameter S in equation (20) is calculated as S' = 0.75. We 
show simulations results, where the charge number per 
particle is Z = 30, 50, 60, 80, 100, and 160. If it is 100 
and the radius R is 25nm, the surface electric field Eg 
becomes 100 V/^m. We also set Ep = Eq. In the case of 
one particle, the interior dielectric constant Ep does not 
affect the exterior electric potential and is irrelevant. In 
the case of two particles, we also performed simulation 
with Ep — 2eo in the examples in figures 6 and 8, but no 
marked difference was found. 

In Subsections 4.1 and 4.2, we will neglect the short- 
range anchoring interaction and set w = 0, focusing on 
the electric field effect on the director field. In Subsec- 
tion 4.3, we will include the short-range anchoring in- 
teraction around a charged particle. In our Landau-de 
Gennes scheme, the orientation order is almost uniaxial 
outside the defect cores both for £i > and E\ < 0. 
Thus we will display the director n around the parti- 
cles. Tube-like surfaces in Figures 4-10 will be those where 
fgd^/T = {dVkQij)^/2 = 0.2. This threshold is so high 
such that the resultant tubes enclose defects. 

In addition, we confirmed that the eigenvalues of Xij 
in equation (22) were kept to be positive everywhere in 
the system. For example, in the uniaxial state with Ei = 
I.Seq and S — 0.75, the eigenvalues of Xij: are given by 
X|| = 0.27 and x± = 0.024. 

4.1 A single particle in nematic liquid 

We fisrt consider a single charged particle for the two 
cases, £i > and Ei < 0. Its charge number Z is in the 
range [60,160]. The orientation tensor Qij is independent 
of the sign of Z. 
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Fig. 4. (color online) Orientational fields around a charged 
particle of (a) Z = 60, (b) Z = 100 and (c) Z = 160 in a 
nematic solvent with £i = — l.Seo- 




Fig. 5. (color online) Top and side views of the director field 
n = (nx,ny,nz) around ansae, corresponding to Z = 160 in 
the panel (c) in Figure 4. 
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small system, the Saturn ring is confined within the box. 
However, if the system size is larger, the defect should be 
more extended, since its radius is predicted to be of order 
£ in equation (24) [23] . 

Figure 4 displays the orientation field around a particle 
with R = 12. Sd for ei = — LSeq- The other parameters are 
the same as in Figure 3. For ei < 0, the director tends to 
be along the particle surface, analogously to the case of a 
neutral particle with planar anchoring realized for w < 0. 
For not large Z in (a), the director is distorted around 
the particle without defects. Slightly above the threshold 
in (b), defects are formed at the two poles of the particle. 
For a neutral particle, a similar defect structure is called 
"boojum" [TirfllS^. For larger Z in (c), two "ansa"-shaped 
defects emerge with their ends on the particle surface, 
as a novel defect structure. Here a boojum-like structure 
in (b) grows into a curved disclination line of topological 
strength s ~ —1/2. The director is perpendicular to the 
plane formed by each ansa. In Figure 5 we show the top 
and side views of the director field around the ansae at 
Z = 160 in Figure 4. Here the axial symmetry is broken, so 
the previous simulation did not detect this structure [33] 
(where an axially symmetric, biaxial defect was instead 
detected) . 

For ei > it was shown [23j that a Saturn-ring ap- 
pears discontinuously with increasing €(oc Z) in equation 
(24) . Also in the case of a neutral particle [3] , its appear- 
ance is discontinuous with increasing wR. In Figure 6, we 
show the normalized gradient free energy J-g/T versus Z 
for Ei/sq — ±1.8, since J-g in equation (4) is sensitive 
to the defect formation. The arrows indicate the point 
of the defect formation on the curves. Remarkably, for 
£i = I.Seo > 0, J^g jumps at Z = 136 with increasing Z 
and at Z = 80 with decreasing Z, where £/R = 5 and 9, 
respectively, using equation (24). This hysteretic behavior 
demonstrates that the system is bistable with and with- 
out a Saturn ring (in the range of 80 < Z < 136 in the 
present example). On the other hand, for ei < 0, J^g in- 
creases smoothly as Z increases. This is because the ansa 
defects gradually protrude from the particle surface into 
the liquid crystal. 



Fig. 6. Normalized gradient free energy J^g vs charge number 
Z of a particle. Arrows indicate the point of defect formation. 
A jump appears for ei > 0, while there is no jump for ei < 0. 

Figure 3 displays the director field n = (nx,ny,nz) 
around a single particle with Z — 60, 100, and 160. Here 
we set El — 1.8£oi and R ~ 12.5d. The liquid crystal is 
deeply in the nematic phase. At the particle surface one 
of the perpendicular alignment is selected, which is anal- 
ogous to the case of a neutral particle in the homeotropic 
anchoring condition realized for w > 0. The system is 
axisymmetric, as assumed in our previous simulation p3 . 
For small Z in (a), no defect is formed, while the orienta- 
tion field is largely distorted. For large Z in (b) and (c), 
a Saturn-ring disclination line of the topological charge 
s = —1/2 appears near the equator of the particle. In our 



4.2 A pair of charged particles 

We first place a pair of positively and negatively charged 
particles with R = 6.25, which form a dipole. Their dis- 
tance is fixed at \Ri — i?2| = 2i?. In Figure 7, we show 
snapshots of the director and the defect structure around 
the two particles. Here Si = LSeq and Zi = —Z^ = 50 
in the upper plates (a) and (b), while E\ = — 1.8eo and 
Z\ = —Zi = 100 in the lower plates (c) and (d). The par- 
ticles are aligned in the parallel direction (left) and in one 
of the perpendicular directions (right) with respect to the 
background director direction (along the z axis for £i > 
and along the x axis for Ei < 0). We can see Saturn rings 
in (a) and (b), while there are four ansae in (c) and two 
ansae in (d) . In the lower panel of Figure 8, we show the 
sequence of this topological change of the defect structure 
with varying the angle between the background director 
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Fig. 7. (color online) Director and defect structure around 
oppositely charged particles in the parallel direction (left) and 
in one of the perpendicular directions (right) with respect to 
the background director direction. Upper plates (a) and (b): 
£i = l.Seo and Zi — —Z2 = 50, where the free energy is lower 
for (a) than for (b). Lower plates (c) and (d): ei — —l.Seo and 
Zi = —Z2 = 100, where the free energy is lower for (d) than 
for (c). 



and the vector connecting the particle centers. In the up- 
per panel of Figure 8, we show the free energy !F — T{9) 
measured from its minimum J-^(0min) as a function of 9. 
The angle ^min at the minimum is for ei > and 7r/2 
for £1 < 0. 

We next place identically charged particles separated 
by 2R. In Figure 9, we display the defect structures around 
two positively charged particle with Z\ = Z^. Remarkably, 
the topology of the defects around a pair is the same as 
that of a single particle. That is, we find only one discli- 
nation loop for ei > and two ansa defects for £\ < 0. 
Notice that a pair may be regarded as a non-spherical 
particle [10] with charge 2Ze. Figure 10 displays the free 
energy = J-{0) measured from its minimum as a func- 
tion of the angle 9. The angle ^min at the minimum is 7r/2 
for £1 > and for Si < 0. 



Fig. 8. (color online) Upper panel: Free energy difference 
AJ^ = J-{6) — J-{Omin) for oppositely charged particles as a 
function of the angle 9 between the background director and 
the vector connecting the two particles. Here Zi = —Z2 = 50 
and 6niin ~ for the curve of £1 > 0, and Zi = —Z2 ~ 100 and 
Omin — n/2 for the curve of ei < 0. Lower panel: Topological 
changes of the equilibrium defect structure for ei — —l.Seo for 
fixed 9 — mr/20 {n — 0, ■ ■ ■ , 10), corresponding to the lower 
panels of Figure 7. 

4.3 A charged particle with nonvanishing w 

In this subsection, we discuss the effect of the short-range 
anchoring free energy Ta in equation (6) supposing a sin- 
gle particle. As illustrated so far, the electric field for pos- 
itive and negative £1 serves to induce homeotropic and 
planar alignment, respectively. Therefore, the two anchor- 
ing mechanisms can compete for (i) £1 > and w < and 
for (ii) £1 < and w > 0. 

In Figure 11, we set Si — LSeq > and choose vari- 
ous negative w. In (a), fgd^/T = 0.03 on the surfaces (in 
green). This threshold is small and the right two snapshots 
do not involve defects. We can see that the region having 
large fg moves from the vicinity of the Saturn ring to up- 
per and lower surface parts of the particle. The Saturn 
ring remains nonvanishing for small w, but the director 
field around the equator tends to be tangential to the sur- 
face (parallel to the background director direction along 
the z axis) and the Saturn ring disappears with increas- 
ing \w\. The director field changes steeply near the sur- 
face away from the equator for large \w\. This changeover 
occurs discontinuously with sudden disappearance of the 
Saturn ring at w = Wc, where Wc = ^iATdr^ in the 
present case. 

In Figure 12, we set £1 — — 1.8£o < and choose var- 
ious positive w. In (a), the ansa defects shrink into two 
point defects and disappear with increasing w. The top 
and side views of the director are shown for w ~ 2Tdr^ 
in (b) and for w = 5Td~^ in (c) around the particle. For 
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Fig. 9. (color online) Director and defect structure around 
identically charged particles in the parallel direction (left) and 
in one of the perpendicular directions (right). Upper plates (a) 
and (b): ei = l.Seo and Z\ = Z2 = 30, where the free energy is 
lower for (b) than for (a). Lower plates (c) and (d): ei = — 1.8eo 
and Z\ = Z-2 = 80, where the free energy is lower for (c) than 
for (d). The defect topology is the same as in the single particle 
cases in Figures 2 and 3. 
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Fig. 10. Free energy difference AT = T(e) - for 
identically charged particles as a function of Q. Here Z\ = Z-i, = 
30 and Smin = 7r/2 for the curve of £1 > 0, while Z\ = Z2 = 80 
and ^min = for the curve of e\ < 0. 
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Fig. 11. (color online) Results at ei = 1.8eo for negative w. 
Surface of fgd'^/T = 0.03 for wd^ /T = -2, -3, -4 and -5 in 
(a). Director n in the xy plane {z — 32d) (left) and in the 
xz plane (y = 32d) (right), where w = ~Td~^ in (b) and 
w = —5Td~^ in (c). The charge and radius of the particle are 
Z = 160 and R = 12. 5d. 



large w, the defect structure becomes axisymmetric with- 
out defects and the regions of large fg covers the particle 
surface. This crossover is continuous with increasing w. 

In the above examples, the short-range anchoring is ef- 
fective close to the surface for sufficiently large \w\, while 
the electric-field anchoring is dominant far from the sur- 
face. A similar problem is encountered in the Fredericks 
transition in magnetic field as the strength of the surface 
anchoring is varied [25,33 . The crossover from weak to 
strong short-range anchoring occurs for 

\w\ > L/Cc - T/d^,, (37) 

where is the thickness of this transition layer. For |ei| < 
Eo this length is determined by 

= E^Jeil/SnLS, (38) 

where Eg — Ze/eoR"^ is the surface electric field. This 
estimation is obtained from hij = in equation (17). On 
the right hand sides of equation (16), the gradient term 
becomes —LS'S/'^ip in the nematic phase, where (p is the 
angle of the director with respect to the surface normal. 
The balance of this term with the last electrostatic term 
(~ eiE-^/8n) at the surface yields equation (38). For our 
parameters chosen in Figures 11 and 12, equation (38) 
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Fig. 12. (color online) Results at ei = — l.Seo for positive w. 
Surface of fgd^/T = 0.03 for wd^ /T = 2,2,, A and 5 in (a). 
Director n in the xy plane [z — 32d) (left) and in the xz plane 
{y = 32d) (right), where w = 2Td'^ in (b) and w = 5Td'^ in 
(c). Top and side views in (c) are indistinguishable. The charge 
and radius of the particle are Z — 160 and R = 12. 5d. 

gives = 0.91fi and the right hand side of equation (37) 
becomes l.lTc?"^, which are consistnt with our numerical 
results. 



5 Summary and Remarks 

We have performed three dimensional simulations in the 
presence of charged particles in nematic liquid crystals. 
We first give a summary. 

(i) The director tends to be parallel (perpendicular) to 
the electric field for positive (negative) ei. In Figure 3, a 
Saturn-ring defect is formed as £ in equation (24) much 
exceeds the particle radius R. In Figures 4 and 5, we have 
found novel ansa defects without axial symmetry in a ne- 
matic solvent with Si < 0. In our previous simulation [23 , 
a boojum-like defect was derived for ei < 0, since it was 
based on the assumption of axial symmetry. In Figure 6, 
the formation of a Saturn ring due to electric field is first- 
order, while that of ansa defects is continuous. 

(ii) We have also examined the director in the presence 
of two charged particles in nematic liquid crystals. Re- 
sults for Zi = —Zi are in Figures 7 and 8, while those 
for Z\ = Zi are in Figures 9 and 10. We have found that 
oppositely charged particle pairs are likely to be aligned in 
the parallel direction for ei > and in the perpendicular 
plane for ei < with respect to the background director 



direction. We conjecture that polar molecules composed 
of oppositely charged particles can be aligned in nematic 
liquid crystals even on microscopic scales. On the other 
hand. Figure 10 shows that the preferred alignment direc- 
tions are exchanged for identically charged particles, 
(iii) We have examined competition of the charge-induced 
anchoring and the short-range anchoring in Figures 1 1 and 
12. These two anchoring mechanisms can compete when 
£i and w have different signs. Under the condition (37), 
the short-range anchoring can be effective near the surface 
with distance shorter than in equation (38). 

We supplement the discussion in Subsection 2.3. For 
microscopic particles (ions), observation of nanoscale de- 
fects should be difficult, but there might be some indica- 
tion of the defect formation in the behavior of the electric 
conductivity [24,,25j. For colloidal particles, the condition 
R > Rc can be satisfied only when the ionization on the 
surface occurs to a sufficient level in a hquid crystal. We 
may also suspend a macroscopic particle in a liquid crys- 
tal. We mention an experiment [30 , in which an electric 
field was applied to nematics containing silicone oil par- 
ticles to produce field-dependent defects. We may even 
propose to suspend metallic particles or water droplets 
containing salt in a liquid crystal, where a surface charge 
appears in an applied electric field. Recently, electric field 
was applied to two-dimensional colloidal crystals in ne- 
matic solvent where the lattice spacing changes up 
to 20% in one direction in response to the applied field. 

Further remarks arc as follows. 

(1) The competition of the short-range and charge anchor- 
ing mechanisms should be investigated furthermore, since 
our examples of a single particle are very fragmentary. 
The interaction among charged particles in liquid crystal 
solvent should be much complicated in such situations. 

(2) The liquid crystal order S increases with increasing \A\ 
in the nematic phase and its discontinuity at the transi- 
tion decreases with decreasing B in equation (2) . For small 
B (for weakly first order phase transition), therefore, the 
defect formation takes place considerably far below the 
nematic-isotropic transition. The ion mobility in nematics 
[24,25j might decrease discontinuously at the Saturn-ring 
formation with lowering the temperature. 

(3) Light scattering should be sensitive to doped ions 
in nematics, where even a small amount of ions should 
strongly distort the nematic order. This is analogous to 
the role of microemulsions in nematics |171I18| . 

(4) As discussed below equation (26), we should examine 
how ionic surfactant molecules can be attached to sur- 
faces of colloids and microemulsions in liquid crystal sol- 
vents [34j . It is worth noting that an ionic surfactant was 
attached to microemulsion surfaces in the previous exper- 
iment [T7]. 

(5) Intriguing also are effects of salt at weakly first order 
nematic-isotropic phase transition and the ion distribu- 
tion around isotropic-nematic interfaces. Such theoretical 
studies were already reported for electrolytes with binary 
mixture solvents [251155] . 
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